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Overview 
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e Motivation 


-Physics-Based Impact Risk Model 
-Airourst & Overpressure Threat 
-Surface Impact Threat 
eFrequency of Interactions 


eResults 
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Credit: anosmicovni 
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Credit: Aleksandr Ivanov 
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Asteroid Threat Assessment Project ra 
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Warning Time (Years 
Credit: Tim Warchocki 
Risk-Informed 
Decision Support 


¢ Mitigation Planning 
¢ Defense Strategies 
¢ Response Decisions 
¢ Policy Development 


Surface Impact Effects 
Modeling 

¢ Ground Damage 

¢ Tsunami Propagation 
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Impact Risk Assessment 
Tools 


Physics-Based Impact 
Risk Modeling 

-* Quantitative Risk Metrics 

¢ Sensitivity to Uncertainty 
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Energy Deposition Modeling 


¢ Analytic model of asteroid entry/breakup to estimate energy deposited in the atmosphere 
¢ Combines progressive breakup of independent fragments and “pancaking” debris clouds. 
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Entry flight: integrates meteor equations of motion and ablation 
dm/dt = -0.5p,;,v°Ao 
dv/dt = p,;,v?AC p/m — gsin@ 
d0/dt = (V/(R;,+h) — g/v)cos@ 
dh/dt = vsin@ 





Fragmentation when pressure > strength 


Pai? > Strength 
n fragments —— debris cloud Each break yields: 


see ¢ Multiple independent, identical fragments (baseline 2) 


Se a ¢ Debris cloud of specified mass fraction (baseline 50%) 


Vis Fragment strengths increase with decreased size 
‘l S, = S,(m,/m,)* 
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Clouds broaden and slow under common bow shock 
= 1/2 
Vdispersion i V ond DP airA’P cloud) 


Energy Deposition, 200m End-to-End Test Case 





Energy deposition computed as change in total KE 
eres | of all fragments/clouds as a function of altitude. 


Airburst at altitude of peak energy deposition. 
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Blast Propagation 
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¢ Blast propagation through the Cart3D asteroid entry modeling 
atmosphere, Simulated using Cart3D at eee 
¢ Unstructured, Cartesian, Eulerian ‘ | 
solver “— 
-lsothermal atmosphere 40 
Energy added to flow in multiple ways 36 
— Static spherical charge = 
— Line source 
¢ Time-dependent pressures tracked 15 
on the ground 10 + Mach contours time = 35 sec 


Peak pressure magnitude 
¢ Arrival times 
¢ Computed results compared with 
Popova et. al (2013) 
-~20m diameter 
°~20 km/s entry velocity 
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Blast Propagation 


ra 
= = 
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Local Overpressure 


ed 





d(Lp)/Lp 
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Land Impact Simulation 


¢ Preliminary test: 040m Iron meteor, 15 km/s vertical impact 


Preliminary Test Simulation. 
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“reser Qcean Impact Simulation 


Preliminary Test Simulation. 


Se 






- 2.1 Million years ago in South-East 
Pacific Ocean 


e Only Known impact into deep Ocean 
Basin 


« Evidence of mega-tsunami debris 
on coasts of Chile and Antarctica as 
well as drill cores from 
Bellingshausen Sea. 


« Current estimates @750m rock at 12 


km/s vertical (10 GT) or 18 km/s at 
45° 


« 9000m deep ocean, 250m 
sediment, basalt crust. 
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Asteroid Characterization ©@---------------9-999999---->- > Input Parameter Distributions 
. . ; o | PHA Measurements 
“,..:.? a a a gi ANT ¢ H-magnitude 
ee? a, “we (eee pansHSeeeee- Monte Carlo Sampling © — - Albedo 
"e fet Ys rat 3 ft ¢ Orbital trajectory 
se yy ¢ Asteroid class 
. Initial cg M4 ¢ Composition 
y Conditions NG... | Impact Parameters 
: Sy ball ¢ Diameter 
a ¢ Density 
« °* Flight Integration > Strength 
(meteor equations of > Luminous efficiency 
motion, ablation) , Fragment-Cloud Model - Velocity 
_ \ (breakup and energy deposition) * Entry angle 


¢ Azimuth angle 


~~ Airburst Altitude ~~ ~~~ . ———— 


(peak energy deposition) 


Local Land Impact Casualties 
(Gridded population within largest damage area) 
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2.5’ (4.6 km) 













































































2.5’ (0.6-4.6 km) 


Global Effects Casualties 
(Percentage world population 
killed by climatic effects) 
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Bolide Events 1994 — 2013 


Small Asteroids that Disintegrated in Earth’s Atmosphere 
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Day (255) _ Night (301) 1 10 100 ania 10000 100000 1000000 


This diagram maps the data gathered from 1994-2013 on small asteroids impacting Earth’s atmosphere and disintegrating to create very 
bright meteors, technically called “bolides” and commonly referred to as “fireballs”. Sizes of orange dots (daytime impacts) and blue 
dots (nighttime impacts) are proportional to the optical radiated energy of impacts measured in billions of Joules (GJ) of energy, and show 


the location of impacts from objects about | meter (3 feet) to almost 20 meters (60 feet) in size. 
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A Near-Earth Asteroid Census se cacy case 


i : aicted total | te " : } @ 
Each image represents 100 objects 
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Credit: NASA 


21 May 2016 ISDC2016, D. Mathias 





Page 13 


Earth Impact Frequency 
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9 Cumulative Impact oe Inputs | 


v. 2015 Inp 
ee 2016 inputs 


Impact Frequency (per year) 





10° 10° 


Min Diameter (m) 


10 10 


Source: Alan Harris 
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Z>. _E¢: Estimated Casualties/ 





Year 





¢ Average potential casualties from impacts within a given size range, 
multiplied by the expected impact frequency/year of that size range. 


E. for each given size range Cumulative E, for impacts 
(changes with bin size) up to given size 
Ec by Size Bin Cumulative Ec 
ae Stony Asteroids, Local Damage 15 Stony Asteroids, Local Damage 
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¢ Expected values based on averages misrepresent likely outcomes 


Conditional Casualty Probability Histogram, Stony Asteroid Impacts 


Conditional Probability 
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act of Population Distribution 
and Global Effects 


¢ Impact location dominates until global effects phase in. 
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Global 

~1500m: Local Pop. 
Avg. Pop. 

Global -— 


~900m: Local Pop. + 
Avg. Pop. 
Global = a. 
~150m: Local Pop. -——${——$—$ $< _} 
Avg. Pop. + 
Global 
~f/0m: Local Pop. 
Avg. Pop. 
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Questions? 
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BACKUP 
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¢ Annual probability of an asteroid up to a given size impacting Earth 
and atfecting at least a given population threshold or greater. 


Absolute Damage Exceedance by Size Threshold 
Bottke Cases 4/20/16, Harris Total Freqs 4/24/16 
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° e.g., contour color plotted at 100m 
diameter and 10° population gives 
the cumulative probability of an 
asteroid 100m or smaller affecting at 
least 1,000 people. 


¢ Probabilities are shown on log-scale 
labeled by order of magnitude, e.g., - 
6 = one-in-a-million (10°) chance per 
year. 


¢ Cumulative probabilities are used to 
minimize result dependence on width 
of diameter and population bins. 


° Affected population represents 100% 
of the population within a 4-psi blast 
damage radius. 





